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Abstract 

We study the influence of converging cooling flow bulk motions on the Sunyaev- 
Zel'dovich (SZ) effect. To that purpose we derive a modified Kompaneets equation 
which takes into account the contribution of the accelerated electron media of the 
cooling flow inside the cluster frame. The additional term is different from the usual 
kinematic SZ-effect, which depends linearly on the velocity, whereas the contribution 
described here is quadratic in the macroscopic electron fluid velocity, as measured in 
the cluster frame. For clusters with a large cooling flow mass deposition rate and/or 
a small central electron density, it turns out that this effect becomes relevant. 
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1 Introduction 



The SZ-effect is becoming more and more an important astrophysical tool 
thanks to the rapid progress of the observational techniques, which allow in- 
creasingly precise measurements. It has thus become relevant to study further 
corrections to it, such as relativistic effects ( [Rephaeli] , |1995| ), the shape of 



the galaxy cluster and its finite extension or a polytropic temperature profile 
(see e.g. [Fuy et a"T| ( |2000[ )), the presence of cooling flows ( |Schlickeiser| , |1991| ; 



[Majumdar and Nath| , |2000|) , corrections induced by halo rotation flCooray and 



Chen| , [200 1|) or even by Brillouin scattering ( (Sandoval- Villalbazo and Maarterii 



2002| ) and the influence of early galactic winds ( [Majumdar et al\ [200 1| ). These 
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additional effects are of different relevance and often depend on the specific 
values of the parameters which describe a given cluster. Generally, they range 
from the percent level up to even 20-30% and accordingly in the subsequent 
determination of the Hubble constant. Taking into account such corrections is 
on one side relevant when determining the Hubble constant via the SZ-effect, 
and on the other hand they could be interesting as a tool to study the detailed 
structure of the cluster itself. 



Following these lines we study here another effect, which has not yet been 
taken into account: The possible influence on the SZ-effect of the cooling flow 
bulk motion of the electron media inside the cluster frame. Indeed, one ex- 
pects such a motion to be present in cooling flows, for which in some clusters 
there is evidence in the central regions ( |Allen et al. , 1993 ; Fabian et al. , 1991 ; 
[Fabian| , |1994| , |1997|) . Recently, |Dupke and Bregmanj ( |2002|) investigated the 
prospects to detect gas bulk velocities in clusters of galaxies through the kine- 
matic SZ-effect, which depends linearly on the velocity component along the 
line of sight. Bulk motions along a given line of sight will contribute to the 
kinematic SZ term as long as their averaged velocity, in the observer frame, 
does not vanish. In the case of cooling flow bulk motion the averaged velocity, 
in the cluster frame, along a given line of sight vanishes in good approximation, 
since we assume spherical symmetric infall. Thus, the cooling flow bulk mo- 
tion will not contribute as such to the kinematic SZ-effect. Indeed, the effect 
we consider depends quadratically on the velocity and clearly, the averaged 
quadratic velocity does not vanish along a line of sight in the cluster frame. 
As we will see, the considered effect usually turns out to be small, since the 
cooling flow bulk motion velocities are rather small, unless in the very central 
regions of a cooling flow clusters, where the cooling flow approaches a sonic 
radius and changes from the subsonic to the supersonic regime. Nonetheless, 
in some favourable cases it might be of the order of some percent of the ther- 
mal SZ-effect. 



The aim of this paper is to examine the influence on the thermal SZ-effect of 
the moving electron media inside a cooling flow region of a galaxy cluster. As it 
will be seen in more detail later, this effect has clearly to be distinguished from 
the kinematic SZ-effect, which takes into account the motion of the whole clus- 
ter or a fraction of it along a given line of sight. In fact, the standard thermal 
SZ-effect describes the frequency dependent intensity change of the Cosmic 
Microwave Background (CMB) photons by inverse Compton scattering off 
the hot intracluster plasma (electrons), where the electrons are supposed to 
be static scatterers. 



The paper is organised as follows: In section 2 we briefly outline the dynamics 
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of a homogeneous steady-state cooling flow model in order to get the velocity 
profile of the bulk motion as well as the corresponding electron density distri- 
bution. In section 3 we derive the modification to the standard Kompaneets 
equation due to the inclusion of bulk motion of the scatterers inside the cluster 
frame and in section 4 we present and discuss our results. In section 5 we give 
a summary and an outlook. 



2 Velocity profile for a homogeneous cooling flow model 



As an example to get a velocity profile which describes the bulk motion, we 
consider a homogeneous steady-state cooling flow model where the mass de- 
position rate m is constant, negative and enters as a parameter in the model. 
No mass drops out of the flow. We neglect the possible influence of magnetic 
fields, rotation and viscosity. In this context, the cluster is expected to be in 
a relaxed state, so that hydrostatic equilibrium allows us to use an isothermal 
/5-model ( Sarazin , 1988|) . For spherical symmetry, the cooling flow can thus 
be described by a set of Euler equations. Mass, momentum and energy con- 
servation read ( [Mathews and Bregman| , |1978| ; |Fabian et al.| , |1984| ; |White ancl 
I987a|jfe|; EarazmL |198% 



B 



arazm 



47rr 2 p CF v = rh 

dr pcf or T 
dE _ Pv dp CF _ a 



where Pcf(j) and P{r) are the gas density and pressure, respectively, in the 
cooling flow and v(r) < is the velocity of the inward directed cooling flow. 
M.{r) is the gravitating cluster mass inside the radius r, discussed earlier in 
the literature by the above authors, with the dark matter {DM) density profile 
as follows: 

Pdm{t) = ?°, , 2 , (2) 
1 + {r/r c y 



where the central cluster density p = 1-8 ■ 10~ 25 g cm' 3 and the core radius 
r c = 250 kpc describe the profile of the cluster mass. As usual, we assume 
that the cooling flow makes no significant contribution to the cluster mass 
density and, that the gas self-gravity can be neglected. The internal energy is 
E{r) = §0(r), with the temperature parameter 9 which defines the square of 
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the isothermal sound speed c s : 

9(r) := c 2 s (r) 



kT e>CF (r) 
\im v 



(3) 



where \i ~ 0.63 is the mean molecular weight, m p the proton mass and T e>OF the 
electron temperature in the cooling flow. A(6) defines the radiative optically 
thin cooling function for a low density plasma, which is defined so that Ap 



CF 



is the cooling rate per unit volume in the gas. We take an analytical fit which 
has been used in earlier works ( [Sarazin and White] , |1987| ; [Majumdar and Mathj , 
200q ): 



A(9) 



10 22 erg cm 3 s 1 



+ 0.313 • T 08 -exp 



4.7 • exp 
T 



T 



3.5 • WK 
i.i 



4.5 



-0.2 



+ 6.42 • T — ■ exp 
+ 0.000439 • T 35 . 



3.0 • 10 6 K 

rp x 4.0' 

2.1 ■ lO 7 ^ 



(4) 



Eliminating the density p CF with the mass conservation equation, the Euler 
system in (1) leads then to a system of two ordinary coupled differential equa- 
tions for the cooling flow velocity v(r) and the square of the isothermal sound 
speed 6{r) ([Mathews and Bregman| , |1978| ): 



§ = v ( 3GM 



lOr0 + f^)/(r 2 (50-3^ 



f = 2 (e{2v 2 r - GM) - {v 2 - 



l m A(<>)\ 



r 2 (56 - 3v 2 ) 



(5) 



To integrate this system we have to take into account that both equations 
have singularities at the sonic radius r s , where 58 (r s ) = 3v 2 (r s ) and where 
the cooling flow undergoes a transsonic transition. The cooling flow region is 
limited by the cooling radius r coo i w 100 — 150 kpc. 

We stress that our goal is not to develop a sophisticated cooling flow model, but 
to get an idea of how the cooling flow bulk motion contributes to the SZ-effect. 
We, therefore, do not attempt to find solutions to the system of Eqs.(5) for the 
innermost supersonic region, r < r s , where we expect shocks to be important. 
We then avoid the search of critical values^ at r s which would have to be 
matched to hydrostatic equilibrium at r cooZ . We thus start our integration from 

1 A transition from subsonic to supersonic occurs, if the numerators and denomina- 
tors in the system of Eqs.(5) vanish at r s . Since the expressions are indeterminate, 
they have to be replaced by nonsingular expressions by the Bernoulli-de l'Hopital's 
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r CO oi towards r s and we stop when the Mach number M = v/c s is close to unity, 
M ~ 0.9. Reasonable initial values at r coo i are: v(r coo i) « v t ~ 10 — 20 km/s, 
where is the turbulent velocity and 9{r coo i) such that t coo \ = f pc ^, A < 
tnubbie- in Fig.l we have plotted the velocity |f(r)| and the isothermal sound 
speed c s (r) as a function of the radius r. This figure illustrates clearly that 
the electrons are strongly accelerated in the cluster center. (For an extended 
discussion: see [Koch| ( |1999| ).) Fig. 2 shows the corresponding electron density 
profile ncF{ r ) i n the cooling flow region. 



3 SZ contribution due to cooling flow bulk motion 



It is well known that the Kompaneets equation ( [Kompaneets] , |1957[ ), result- 
ing from a Boltzmann transport equation, describes the spectral shift of the 
CMB photons due to the inverse Compton process on the hot electrons of the 
intracluster medium: the SZ-effect - for a review, see pirkinshawj ( |1999| ). 
The frequency dependent intensity change of the CMB photon field AIk(x) 
after integration along the line of sight over the cluster (cl) dimension can be 
expressed as follows ( |5unyaev and Zerdovich| , |1972| ; [Rephaeli| , |1995|) : 



A/* (re) 



i g(x 




a T n d dl c 



(6) 



where x = ^ is the dimensionless frequency with T the CMB temperature 

and i = 2 n^X ■ The integral is the Comptonization parameter y K describ- 
ing the cluster properties with T ejC z,m e the electron cluster temperature and 
mass, respectively. n c \ is the electron number density in the cluster and <jt 
the Thomson cross section, (k, h, c are the Boltzmann constant, the Planck 
constant and the speed of light, respectively.) The function g(x) defines the 
spectral shape of the thermal SZ-effect where, for the plasma in clusters, we 
have Tl d > T: 



x 4 expx /x(exp x + 1) 
(expx — l) 2 \ expx — 1 



(7) 



It is crucial to notice that the standard Kompaneets equation describes a static 
scatterer, assuming that in the average the electrons are macroscopically at 
rest. This is no longer true for the electrons in an accelerated cooling flow, as 
it is described by the system of Eqs.(5). The Kompaneets equation has thus to 

rule, which will give initial values for Eq.(5) at r s if one assumes a temperature at 
the sonic radius. 
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be modified in such a way, that the (macroscopic) bulk velocity of the moving 
electron media is explicitly taken into account. 

[Psaltis and Lamb] ( |1997| ) gave a very detailed analysis of Compton scattering 
by static and moving media. They made a careful distinction between the 
electron rest frame, the fluid frame (comoving with the fluid) and the system 
frame. Starting from the Boltzmann equation in the system frame, introduc- 
ing the proper Lorentz transformations, expanding to the appropriate orders 
and assuming that the velocity distribution in the fluid frame is a relativistic 
Maxwellian, they end up with the zeroth moment of the radiative transfer 
equation with emission and absorption included. 

Under the condition that the radiation field (CMB) is isotropic in the sys- 
tem frame, and introducing the macroscopic electron bulk velocity v(r), the 
extended Kompaneets equation becomes (Psaltis and Lambj [19971) : 



1 dl 



n e aTC dt m e c 2 de 




a 91 d \ n 



where we neglected absorption, emission and induced scattering, e is the pho- 
ton energy and I the corresponding intensity. 

For CMB photons we have e <C kT SyC i + |m e t> 2 . Thus, the above Eq.(8) be- 
comes: 

1 dl _ (kT efil | v 2 \ 
n e arc dt \m e c 2 3c 2 J 
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de 2 



(9) 



Eq. (9) shows that if the radiation field is isotropic in the system frame, Comp- 
tonization of the bulk motion of the electrons inside the cluster is described 
entirely by second order terms in v, since all first order terms in v vanish 
identically. The effect is clearly different from the kinematic SZ-effect, where 
the cluster as a whole moves through the CMB radiation. This equation re- 
duces, of course, to the standard Kompaneets equation for v = 0. Thus, the 
important point is that the effect of a non-zero velocity field gives an additive 
contribution to the standard thermal SZ-effect. 

This way we can express the bulk motion contribution due to the cooling flow 
(CF) to the SZ-effect as follows: 

di (di\ (di\ , . 



where the first term on the right hand side is given by the standard Kompa- 
neets equation and the second term is due to the moving electron media in 
the cooling flow. 
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Relating the CMB photon intensity field / to the occupation number r\ through 
/ = ioX 3 r], we have for the cooling flow contribution from Eq.(9): 




CF 



d d 2 
-4x—(x 3 r]) + x— -(x 4 ?]) 
ox ox z 



. o T n CF v 2 

*o (11 

c 3 



where n CF is the electron number density in the cooling flow region. Integrating 
over the cluster cooling flow region and assuming a Planckian photon field for 
i], we find: 

1 f v 2 (r) 
AIcf(x) = i g(x)— I n CF (r)a T ^—dl CF , (12) 



c 

CF 



which is the contribution due to the cooling flow bulk motion inside the cluster 
frame to the SZ-effect. From Eq.(9) we see that the spectral shape g(x) of this 
additional term is the same as for the usual thermal SZ-effect. This clearly 
makes it more difficult to distinguish this contribution from the usual SZ-effect. 



4 Results 



In the following we get an estimate of this additional SZ contribution. In order 
to compare bulk motion contribution in the cooling flow region to the classical 
thermal SZ-effect, we introduce the following frequency independent ratio: 

A/cf Vcf no , 
y K 



with y CF = J n cF(r)v 2 {r) a ^ and = j ^e^ aTn j r \ ^ Comptoniza- 

cf 60 d meC 

tion parameter. 

We consider a typical spherical cluster with an isothermal /?- model for the 
electron cluster density n c i. When the line of sight goes through the center of 
a spherical cluster, we get for a large cluster radius r d , r cl 3> r c - see e.g. ( [Puy 
|it~aTl , |2000l) : 



d^¥-?2^ (14) 
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where r c is the core radius, no,d the central electron density of the cluster and 
B the Beta function. This way we find: 



^ ' cool 

a = TT^ "iT^ 1 TT / n CF {r)v 2 {r)dl CF . (15) 



2' 2- 



The integral is limited to the cooling flow region from r coo i to the sonic radius 
r s , as explained in section 2, along the line of sight through the cluster center. 
The integration boundary r s depends on the mass deposition rate and the 
initial value conditions^. 

As an example we take as typical values for clusters of galaxies the numbers 



found for A426, A478, A1795, A2029 and A2142 flWhite, Jones and Formarm 



19971 ). All these clusters have a cooling flow with a mass deposition rate of 
rh ~ —300 M Q /yr. The values for (3 and the central electron cluster density 
no,d are from |Mohr et al.| ( |1999| ). With these parameters we can estimate the 
ratio a between the cooling flow bulk motion contribution and the classical 
thermal SZ-effect, see Eq.(15), which turns out to be of the order 10 -5 . As 
we cut out the most inner part of the cooling flow region, from the center 
to the sonic radius r s , this value must be interpreted as a lower limit to this 
effect. Thus at first glance the cooling flow bulk motion contribution is not 
very relevant for SZ observations and for the subsequent determination of the 
Hubble constant. 

Nevertheless, as a = a(T ejd ,n 0jC i, P,m,r c ,r coo i) depends on several observa- 
tionally inferred parameters which can substantially change, we performed 
also a numerical investigation by letting vary the various parameters entering 
into the determination of a. We find that, for instance: 

• For big mass deposition rates m(^ — 3000 M & /yr) and small cooling and 
core radii (~ 100 — 150 kpc) a increases to 10~ 4 , by assuming that the other 
parameters have still average cluster values such as: T e ci = 7.5 keV, (3 = |, 
n 0>c i = 2.2 • 10- 2 cm- 3 . 

• Whereas we expect (3 ~ | and T £)d ^7.5 keV not to vary substantially, the 
value for the central electron cluster density could easily be smaller by two 
orders of magnitude, thus reaching values as low as nQ tC i ~ 10~ 4 cm~ 3 . If this 
is the case then a can be of the order 10 -3 — 10~ 2 . Notice that ucf, since 
it depends on the total gravitating mass as can be seen from the system of 
Eqs.(l), will not change much by varying no )C /. 

• When taking into account the finite extension of a cluster or an aspherical 
gas distribution it turns out that y K gets smaller by as much as 30% in some 
cases, see [Puy et aT| ( |2000| ). Accordingly, a increases by ~ 40%. 



The sonic radius r s is typically between 1 and 15 kpc. 
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Thus, the cooling flow bulk motion contribution to the SZ-effect can reach for 
some clusters the percent level and become relevant for a correct analysis of 
the SZ observations. Especially in view of the rapid progress in the observa- 
tional techniques - which will provide in the near future much more accurate 
measurements - one might then be capable of observing also small deviations 
from the standard thermal SZ-effect. This will then enable observers either to 
obtain better determinations of the Hubble constant or more details on the 
state of the cluster. 



5 Summary and Outlook 



We showed how the Kompaneets equation has to be modified to include the 
effect of the bulk motion due to a homogeneous steady-state cooling flow. This 
effect is quadratic in the cooling flow velocity, as it is measured in the cluster 
frame, and thus different from the usual kinematic SZ-effect, which is linear 
in the velocity. 

The contribution strongly depends on the specific dynamics of the cooling 
flow and the cluster properties. Generally, it turns out to be rather small. 
However, for some clusters which have more extreme parameters, though not 
irrealistic, one finds a contribution at the percent level. One might speculate 
that this effect could be relevant when observing young clusters which are not 
yet virialized, still in their formation phase and for which one might expect the 
presence of rather large and extended regions with bulk motions, or perhaps 
in superclusters which are still in the formation phase. Such young clusters 
might be observable in future with the Planck satellite. 

Another interesting situation might arise from shocks in merging processes of 
clusters of galaxies, in situations where the average velocity along a line of 
sight approximately vanishes. If so, the merging process will not be detected 
through the kinematic SZ-effect which is linear in the velocity, but could be 
instead observed through the above discussed effect which is quadratic in the 
velocity. 

Acknowledgments 

We would like to thank S. Majumdar and S. Schindler for useful discussions. 
This work was supported by the Swiss National Science Foundation. Part of 
the work of D.P. has been conducted under the auspices of the D r Tomalla 
Foundation. 



9 



References 



Allen, S.W., Fabian, A.C., Johnstone, R.M., White, D.A., Daines, S.J., Edge, 

A.C., Stewart, G,C, 1993. MNRAS 262, 901 
Birkinshaw, M., 1999. Phys. Rep. 310, 97 
Cooray, A., Chen, X., 2002. ApJ 573, 43 
Dupke, R.A., Bregman, J.N., 2002. ApJ 575, 634 
Fabian, A.C., Nulsen, P., Canizares, C, 1984. Nature 310, 733 
Fabian, A.C., Nulsen, P., Canizares, C, 1991. AA Rev. 2, 191 
Fabian, A.C., 1994. Ann. Rev. Astrophys. 32, 277 

Fabian, A.C., 1997. in: Galactic and Cluster Cooling Flows , ASP conference 
Series 115 

Koch, P.M., 1999. Diploma Thesis ETH-Ziirich, unpublished 

Kompaneets, A.S., 1957. Sov. Phys. JETP 4, 730 

Majumdar, S., Nath, B., 2000. ApJ 542, 597 

Majumdar, S., Nath, B., Chiba, M., 2001. MNRAS 324, 537 

Mathews, W.G., Bregman, J.N., 1978. ApJ 224, 308 

Mohr, J.J., Mathiesen, B., Evrard, A.E., 1999. ApJ 517, 627 

Psaltis, D., Lamb, F.K., 1997. ApJ 488, 887 

Puy, D., Grenacher, L., Jetzer, Ph., Signore, M., 2000. AA 363, 415 

Rephaeli, Y., 1995. Annu. Rev. Astron. Astrophys. 33, 541 

Sandoval- Villalbazo, A., Maartens, R., 2001. astro-ph/0105323 

Sarazin, C, 1988. in: X-ray emission from clusters of galaxies, Cam. Uni. Press 

Sarazin, C, White, III R.E., 1987. ApJ 320, 32 

Schlickeiser, R., 1991. AA 248 , L23 

Sunyaev, R.A., Zel'dovich, Ya.B., 1972. Comm. Astrophys. Space Phys. 4, 173 
White, D., Jones, C, Forman, W., 1997. MNRAS 292, 419 
White, III R.E., Sarazin, C, 1987a. ApJ 318, 612 
White, III R.E., Sarazin, C, 1987b. ApJ 318, 629 



10 




Fig. 1. Cooling flow velocity \v(r)\ and isothermal sound speed c s {r) as a function of 
radius. The mass deposition rate is m = — 300 M yr -1 , r coo i = 100 kpc. The initial 
values for the integration are v(r coo i) = 15 km/s and c s {r coo {) = 1600 km/s. 
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Fig. 2. Cooling flow electron number density ncF( r ) as a function of radius. Cooling 
flow parameters as adopted in Fig.l. 
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